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We review the theoretical and experimental situation for long-lived heavy quarks, or bound states 
thereof, arising in simple extensions of the Standard Model. If these particles propagate large dis- 
tances before their decay, they give rise to specific signatures requiring dedicated analysis methods. 
In particular, vector-like quarks with negligible couplings to the three known families could have 
eluded the past experimental searches. While most analyses assume prompt decays at the produc- 
tion vertex, novel heavy quarks might lead to signatures involving displaced vertices, new hadronic 
bound states, or decays happening outside of the detector acceptance. We perform reinterpretations 
of existing searches for short- and long-lived particles, and give suggestions on how to extend their 
reach to long-lived heavy quarks. 



I. INTRODUCTION 

Over the past decades, we have discovered that Nature 
consists of a given number of elementary particles, deeply 
connected with the known fundamental forces driving our 
universe. The recent observation of a new particle resem- 
bling the long-sought Higgs boson at the Large Hadron 
Collider (LHC) now provides us with strong evidence for 
the validity of the Standard Model (SM) [HE] . Yet, while 
it is generally acknowledged that the latter comprises 
three generations of chiral quarks and leptons, various 
fundamental problems do not find their answers within 
this framework. For instance, long-standing issues such 
as the origin of the fermion mass hierarchy or the na- 
ture of the Cabibbo-Kobayashi-Maskawa (CKM) mixing 
matrix hint at the possible need for new physics. While 
the SM is expected to loose its predictive power as the 
experimental searches now reach the TeV scale, models 
suggesting new heavy fermions might offer promising so- 
lutions to these problems in the near future. 

Although a fourth sequential family of quarks is now 
claimed to be disfavoured by the recent results of the 
searches for the Higgs boson [3i , other models predict- 
ing new heavy particles beyond the top quark are still 
consistent with the current experimental measurements. 
In particular, the possibility for vector-like quarks, i.e., 
quarks having their left- and right- handed components 
transforming identically under the electroweak gauge 
group, are a common feature of many scenarios going 
beyond the SM. Grand unified models based on Eq, for 
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example, predict the existence of isoscalar, charge — l/3e 
quarks with small, quite possibly tiny mixing to the well- 
known doublet quarks [1]. Non-chiral quarks also have 
the peculiarity to decouple in the heavy mass limit, lead- 
ing to SM-like signals. Should the mixings of these with 
the light SM fermions be suppressed, the Higgs produc- 
tion rates would not be easily distinguishable from the 
SM expectations [5]. These new quarks could be suffi- 
ciently long-lived to behave like effectively stable parti- 
cles, evading the current searches as they propagate over 
sizeable distances. 

From the collider searches it is clear that if new states 
with masses less than a hundred GeV had existed, they 
would have been observed. On the other hand, the exper- 
imental reach for detecting novel heavy stable particles 
above a few TeV depends on how readily identifiable such 
states are at the LHC. While long-lived particles could 
fly away from the primary interaction point before they 
decay to ordinary particles and lead to displaced vertices, 
they can also hadronise, allowing for a possibly rich spec- 
trum of new exotic and heavy bound states. Such stable 
massive particles are also anticipated at the TeV scale in 
many new physics models, either due to the presence of 
a new conserved quantum number (e.g., R parity in su- 
persymmetry) , or because the decays are suppressed by 
kinematics or couplings (see \§\ and references therein for 
an exhaustive review) . Although the majority of the past 
studies focused on stable states in supersymmetry, the 
possibility for new stable quarks received few attention 
at the LHC. While the CMS and ATLAS experiments 
are now setting strong limits on long-lived gluino, stop 
and stau pair production (see, for instance, [7] and [8]), 
scenarios involving quasi-stable heavy quarks have been 
barely investigated at the LHC. 
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In the following, we review the phenomenology of new, 
long-lived heavy quarks. Assuming a general parametri- 
sation, Section [IT] first examines their production and 
decay modes at the LHC, considering both the chiral 
and vector-like scenarios. The possible signatures for 
displaced vertices and stable massive hadrons are then 
covered within the context of negligible mixings with 
the SM fermions. The experimental aspects of the col- 
lider searches for new long-lived quarks are presented in 
Section |III| We review the situation for prompt decay 
searches, displaced vertices and Heavy Stable Charged 
Particles (hereafter, HSCPs) and highlight some of their 
limitations. Reinterpretations of the existing searches are 
then presented for short- and long-lived particles. Con- 
sidering the possibility to improve the sensitivity of the 
existing analyses to long-lived quarks beyond the TeV 
scale, some alternative search topologies are finally de- 
scribed. Our conclusions are given in Section [IV| 

II. LONG-LIVED QUARKS SIGNATURES 
A. Production 

As we review the possibility for long-lived quarks in 
general, two different scenarios will be distinguished in 
this work. New chiral fermions, for which the left- and 
right-handed chiralities have different charges under the 
electroweak gauge group, gain their masses from the 
electroweak symmetry breaking mechanism. If heavy, a 
sequential fourth generation quark doublet {t' can 
therefore induce large corrections to loop observables. 
Such a scenario is now severely constrained from the 
flavour and electroweak precision data, and from the 
Higgs searches results. Given that chiral fermions couple 
to the Higgs boson with a strength proportional to their 
Yukawa couplings, they do not decouple from its produc- 
tion, as its rate should increase by a factor of about 9 
due to the additional fermion loops occurring in gluon 
fusion. The unobserved enhancement thus now strongly 
disfavours a fourth generation Standard Model [3]. 

Non-chiral quarks, on the other hand, still provide a 
viable extension to the SM and certainly require a dedi- 
cated discussion. Considering the introduction of vector- 
like SU(3)c colour triplets, new parameters have to be 
added, namely, their masses and all associated couplings. 
Such quantities arise naturally from the corresponding 
SU(2)l x U(l)y invariant Yukawa interactions and Dirac 
mass terms in the SM Lagrangian. The light quark cou- 
plings to the Higgs and electroweak bosons are conse- 
quently modified, along with the introduction of new 
Yukawa and gauge couplings between the heavy and the 
SM quarks. 

Still, it is important to emphasise that such new states 
cannot have arbitrary quantum numbers, since they can 
only mix with the SM quarks through a limited num- 
ber of gauge-invariant couplings. Classifying them into 
SU(2)l multiplets, their Yukawa terms only allow for 
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TABLE I. Vector-like multiplets allowed to mix with the SM 
quarks through Yukawa couplings. The electric charge is the 
sum of the third component of isospin T3 and of the hyper- 
charge Y . Adapted from [5]. 



seven distinct possibilities, i.e., the two singlets, the three 
doublets and the two triplets displayed in Table [TJ 

The production of such new heavy quarks, either chi- 
ral or vector-like, is usually assumed to proceed dom- 
inantly at the LHC through gluon fusion, gg — > QQ. 
Nevertheless, depending on the model at hand, elec- 
troweak single production can also provide an alterna- 
tive mechanism as it is not affected by the large phase- 
space suppression from which pair production suffers. In 
particular, new heavy quarks can be produced singly in 
flavour-changing processes via the electroweak interac- 
tion through q^q/ V* -> qtQ where V = W, Z. A 
comparison between the dominant production modes can 
be found in [101 HI] for fourth generation, and in [12] 
for vector-like quarks. Assuming order unity couplings, 
benchmark cross-sections have been evaluated for various 
mass scenarios and reproduced in Table [H] 

As described in [TJJ, the leading contributions to 
vector-like T and B single production arise from ud — ¥ 
W* -> jT, du -> W* -> jB, qd -> Z* -> jB and 
qu — > Z* — > jT, where q denotes a valence quark par- 
ton and j a generic light quark jet. Unlike QCD pair 



channel y/s = 7 TeV y/s = 14 TeV 

m Q = 900 GeV/c 2 m Q = 1800 GeV/c 2 

pp -> W -> jX 1.4 0.36 

pp W jX 0.037 0.0092 

pp^W^jT 0.61 0.16 

pp -> W* -> jT 0.052 0.013 

pp^Z*^jT 0.43 0.11 

pp -> Z* -> jf 0.025 0.0064 

pp -> W -> jB 0.69 0.18 

pp -> W -> jB 0.089 0.022 

pp^Z'^jB 0.18 0.047 

pp -> Z* -> jB 0.034 0.0088 

pp -» W* -> jY 0.29 0.074 

pp^W-fjY 0.12 0.031 



TABLE II. Electroweak single production cross-sections (in 
pb) as a function of the mass tuq of the heavy quark, assuming 
order unity couplings as computed in |12| . In contrast, strong 
pair production pp — > QQ yields 3.61 (0.63) fb for mq — 900 
(1800) GeV/c 2 at = 7 (14) TeV. 
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production, the above processes however scale with the 
Q — q quark coupling squared, and can thus be strongly 
suppressed if the associated mixings are negligibly small. 
On the other hand, Figs. 2 — 4 in [T2] indicate that 
cr(pp — > QQ) falls off faster than the single production 
cross-sections as soon as rag reaches the TeV scale, while 
the electroweak production channels displayed in Table 
[111 take advantage from the possibly large rriq /n%y en- 
hancement originating from the longitudinal polarisation 
of the gauge boson V. Such a behaviour can also be un- 
derstood from the phase-space suppression arising from 
the presence of one or more heavy quarks in the final 
state, as well as from the decreasing parton luminosity 
for large x values at the LHC. Furthermore, given the 
difference in the PDFs of valence and sea quarks in the 
initial states, new heavy quarks are produced singly at a 
much higher rate than antiquarks for increasing rriQ. 

We point out that the relative rates of the jQ and jQ 
channels strongly depend on the valence quark density 
in the initial state. While the u quark contributes more 



than the d quark, the neutral current contributions are 
weaker than the charged current ones. Finally, we under- 
line the absence of the tree-level processes pp — > Z* — > jQ 
for vector-like quarks involving exotic charges, given that 
they can only interact with other states via tree-level 
charged currents [12] . Although most of the above state- 
ments are model-dependent, we conclude that the large 
running energy of the LHC makes single electroweak pro- 
duction a promising discovery channel when considering 
new heavy quarks searches with tuq > 1 TeV/c 2 . 



B. Decay 

In the search of new fermions, the associated decay 
modes require careful attention. The exact partial width 
for a new sequential heavy quark Q decaying on-shell to 
a light quark q through a charged current can be written 

as 



T(Q -> qW) 
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f 2 (a, /3) = [(1- a 2 ) 2 +/3 2 (1 + a 2 ) - 2/3 4 ] y/[l - (a + f3) 2 ] [1 - (a - 
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where kq 9 denotes the generic Q—q quark coupling, equal 
to the CKM mixing matrix element Vq 9 when consider- 
ing a sequential new family of quarks. Interestingly, the 
width |1| holds for both chiral and vector-like quarks, 
given that all heavy-to-light charged current quark de- 
cays occur to be pure V — A processes with identical 
rates [13] ■ If tuq 3> m q , the above partial width can be 
shown to reach the asymptotic form 



r(Q qW) ~ 170.5 MeV x |« Qg | 



(3) 



Considering a new heavy quark Q decaying exclusively 
via charged current 2-body decays, we evaluate in Table 
|n]the corresponding decay lengths, with ct ~ 2 x 10~ 10 
GeV/r(GeV) jxm and \kq q \ = 1. With increasing mass, 
the probability for a new heavy quark Q to decay weakly 
thus becomes more and more important as the individ- 
ual lifetime decreases. However, small couplings to the 
lighter SM quarks strongly affect this statement as the 
total widths are suppressed by |kq<j| 2 . As we will see 
in Section |II D[ this could lead to a very different phe- 
nomenology as the quarks Q, if very heavy, would form 
bound states and possibly decay hadronically. 
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22.91 (83.85) 
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20.90 (11.14) 
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500 


40.94 (27.86) 
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600 


70.81 (54.52) 


2.82 (3.67) 




700 


112.50 (93.06) 


1.78 (2.15) 




800 


167.97 (147.43) 


1.19 (1.38) 




900 


239.18 (213.57) 


83.62 (9.36) xlO 
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1000 


328.12 (299.46) 


6.10 (6.68) xlO" 
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1100 


436.74 (405.05) 


4.58 (4.94) xlO" 
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1200 
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2.22 (2.33) xlO" 


i 


1500 


1107.5 (1063.79) 


1.81 (1.88) xlO" 
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1600 


1344.1 (1297.4) 


1.49 (1.54) xlO" 


i 


1700 


1612.2 (1562.5) 


1.24 (1.28) xlO" 
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1800 


1913.8 (1861.1) 


1.05 (1.07) xlO" 


i 



TABLE III. Partial widths and decay lengths for a new heavy 
quark Q decaying to light quarks via on-shell charged currents 
Q -» qW, with m, = (173.2) GeV/c 2 and \n Qq \ = 1, as 
computed from Eq. |l]). 

Several conclusions can already be drawn at this stage. 
As far as a heavy fourth quark family i s con " 

cerned, we highlight that the ratios of the charged current 
decay rates to the lighter families only depend on the off- 
diagonal mixing elements. Assuming that the extended 
4x4 CKM matrix is unitary, non-zero values for the 
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CKM mixing elements \VQ q \ imply non-unity branching 
ratios for fourth generation quarks in general. If no as- 
sumption is made on the hierarchy of the Q — q coupling 
strengths to weak currents, no exclusive flavour chang- 
ing transitions should then be preferred for new heavy 
quarks. 

Flavour Changing Neutral Currents lead to similar 
conclusions for models involving new vector-like quarks 
interacting with the SM families. Depending on their 
Yukawa couplings, they can indeed undergo tree-level 
transitions to lighter quarks and W, Z and Higgs bosons. 
Considering vector-like quarks Q — T, B with tuq 3> m q , 
the neutral currents Q — > qZ and Q — > qH can occur with 
a rate of the same magnitude as given by ^ , if they are 
allowed to mix with the light SM quarks q. 

For non-chiral quarks decaying to Wq, Zq and Hq, 
long lifetimes can thus be expected in the \nQ q \ — 



limit. Additionally, it is interesting to notice that T(Q — > 
qZ)/T(Q — > qH) ~ 1 holds in the large mass limit 
mq 3> m q for most scenarios, while the charged-to- 
neutral current ratios T(Q — > qW)/T(Q — > q(Z,H)) ap- 
proach either 2 or depending on the vector-like repre- 
sentation [14]. 

If more than one new heavy quark is considered, and 
if they are allowed to decay into each other, the partial 
width (|3| becomes inaccurate due to the possible compe- 
tition with the heavy-to-heavy transitions Q\ — > QzV^, 
where V — W,Z, H , depending on the model at hand. 
Indeed, if the decay rates of the daughter particles are 
significant with respect to the mass of the parent, the de- 
cay of an unstable particle is allowed through threshold 
effects, even if kinematically forbidden. For both chiral 
and vector-like quarks, real and virtual V emission near 
threshold proceed with the rate 



r(Qi -> Q 2 v { * 
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/ 3 (a,/3, 7 ) = 2 



(1 ~^ )2 , [(1 ~ /3) 2 + x (1 + 0) - 2x 2 ] y/\(l,X,P) 
dx 



[(1 - ax) 2 + 7 2 



(5) 



where \(a,b,c) = a 2 + b 2 + c 2 - 2ab - 26c - 2ac g]. 
While the Q —> qV rates essentially depend on the 
Q — q quark couplings, the heavy quark mass splittings 
\m,Q 1 — toq 2 | drive the strength of the heavy-to-heavy 
transitions [TTJ[TS]. If the quarks Qi an d Qi belong to 
the same weak multiplet with a large mass splitting, the 
Qi Q2V decay could be very rapid and dominate the 
other decay modes if there is no \kq 1 q 2 | suppression. On 
the other hand, the lightest partner can only decay into 
a SM family quark, and is thus plausibly long-lived for 
small mixings. Should this be the case, we will see in 
Section that the latter could hadronise, leading to a dif- 
ferent decay phenomenology. This occurs specifically for 
SU(2)l vector-like singlets, which free quark decays can 
only proceed via non-zero mixing into the lighter gener- 
ations. 

As far as doublet and triplet representations are con- 
cerned, new vector-like partners should be very close 
in mass if one assumes that isospin conservation is re- 
spected. Although it may be broken by higher order 
corrections, such a degeneracy has been shown in |16j to 
require a mass difference between 70 and 1 10 MeV for 
non-chiral doublets, allowing for nanoseconds lifetimes. 
If their couplings to the lighter quarks are negligible, the 
partial widths for -X5/3, T2/3, 5— 1/3 and Y-4/3 quarks 
are thus pausibly small. 

Summarising these results, the range of the observable 
couplings to which long-lived particles searches are sensi- 
tive at the LHC can be evaluated. From the approxima- 



tion Q, the partial decay length for a vector-like quark 
reads 



A ~ 2 cm x 



10- 



2 /500 GeV/c 



r)( 



• 2 N 3 



(6) 



which gives an order of magnitude for the required Q — q 
quark couplings to detect a new heavy quark Q. Assum- 
ing mq < 1 TeV/c 2 and an experimental resolution of 
at least 25 yum, we obtain that the heavy quark decays 
are prompt only if their associated couplings are above 
the I0 -7 level. If |kq 9 | < 10 -7 , the corresponding decay 
lengths could be observed over a few tens of centimeters. 
In the limit of such small mixings with the SM fcrmions, 
it is thus natural to motivate searches for displaced ver- 
tices. 



C. Classification of signatures 

We now summarise the possible situations giving rise 
to the possibly observable signatures of long-lived quarks 
at the LHC. 

The classification given in Table |IV| summarises the 
long-lived quark scenarios corresponding to the small 
mixing scenarios treated in this work, namely, KQ q cou- 
plings between 10 -8 and I0 -10 . Furthermore, if mixing 
with all SM quark families is permitted, direct searches 
should aim at being sensitive to all light quarks q in the 
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TABLE IV. Possible decay signatures for new long-lived 
quarks (ibq = 1 TeV/c 2 ), as discussed in the text. 



observed final states. Given that their interactions are al- 
lowed through arbitrary Yukawa couplings, the searches 
for long-lived vector-like quarks should be as inclusive as 
possible in order to cover the full spectrum of possibilities 
at the LHC. Such a program is important to set new con- 
straints on the Q — > qV^ decay modes, independently 
of any assumptions on the Q — q mixing sector. The 
relative importance of the Q (u,c,d,s)V channels is 
most relevant compared to the Q — > (t, b)V charged cur- 
rents, as no exclusive decay mode should be preferred for 
V = W,Z, H. 




io- 10 io~ 8 icr 6 | | Mr* icr 2 l 

\ K Qq\ 



FIG. 1. Lifetime r in seconds (7/3 ~ 1) for a sequential heavy 
quark Q decaying to the lighter families by emitting a real 
W boson, for m Q = 300 GeV/c 2 (in blue) and m Q = 1800 
GeV/c 2 (in red). The four regions of interest refer to the 
conventions given in Tabic |Tv| 

As illustrated in Figure [I] four cases for new long-lived 
quarks can be distinguished depending on their mixing 
parameters with the SM fermions. Firstly, if all couplings 
with the known SM fermions are larger than the 10~ 7 
level, all heavy-to-light decays Q — > qV are prompt, over 
distances smaller than 25 /im. This defines the short- 
lived scenario (i), with no experimental difference com- 
pared to the direct searches. 

The second scenario, (ii), arises for intermediate decay 
lengths ranging between a few microns and centimetric 
distances. If there is at least one light quark q such that 
IO -7 < \nQ q \ < 10 -9 , displaced events could be observed, 
yet with a partial width suppressed by |ftQg| 2 . Conse- 



quently, if all the couplings but one are below the 10 
level, a single exclusive channel would lead to a prompt 
decay signature with a partial width large enough to al- 
low for short-lived heavy quarks. 

A possible exception, though, is the case for new heavy 
multiplets with sizeable mass splittings, as allowed in ex- 
tra generation models and possible extensions [Tl] [T71 - 
119) . If mQ 1 > rriQi + Wy, the heaviest quark Qi can 
be short-lived and decay semi- weakly to Q 2 V<-*\ On the 
other hand, the lightest partner is likely stable if all its 
decay modes suffer severe suppression. If mQ ± ~ itlq 2 , 
all heavy-to-heavy transitions are suppressed and both 
quarks could be long-lived on detector scales. 

For particles with decay lengths larger than the detec- 
tor dimensions lies the narrow, yet unexplored region (Hi) 
for which no strong limits can be set on long-lived heavy 
quarks. Details will be discussed in Section [TlD| Finally, 
if all heavy quarks couplings with the SM fermions are 
below the 10 -10 level, the stable case (iv) becomes a 
relevant scenario, possibly in conjunction with displaced 
events. We emphasise that the scenarios (i), (ii) and (Hi), 
for which prompt decays and displaced vertices could be 
observed with short, intermediate and large lifetimes re- 
spectively, only call for one of the allowed decay modes 
to fulfill the associated conditions. The stable scenario 
(iv), on the other hand, requires all decay lengths and 
partial widths to lie in the ranges given in Table |IV| 

In any case, the observability of the above signatures 
strongly depends on the allowed decay modes. As we 
will see in the next section, if all Q — q quark cou- 
plings to the SM families are smaller than 10~ 2 , such 
new heavy fermions could hadronise. As a result, the an- 
nihilation decays and hadronic transitions between the 
formed bound states would dominate, while the single 
quark decay events might be suppressed. 



D. Signatures for heavy stable particles : heavy 
quarkonia 

Previously, we have discussed how displaced events 
could occur in the single quark decays transitions in the 
case of small Q — q couplings. In this section, we describe 
the possibility for new heavy quarks to bind into hadrons 
if their lifetime is large enough, providing us with an in- 
teresting alternative for new signal searches at the LHC. 

While new coloured particles with masses heavier than 
the top quark are usually assumed to decay as free par- 
ticles, their possible formation into baryons and mesons 
is an important issue to consider in the case of small 
couplings with the SM fermions. 

Should the partial widths Q —> q(W,Z,H) be sup- 
pressed, the quark Q might be stable and hadronise. 
Assuming that the binding force in a fourth generation 
bound state is of Coulombic nature, the seminal condi- 
tion 

m Q < 125 GeV(100 GeV) x \V Qq \- 2/3 (7) 



G 



was derived in 20 for QQ (Qq) formation. According to 
Eq. Q, 

new heavy quarks form hadronic states if their 
mixing with known quarks is sufficiently small. Mixings 
roughly smaller than 0.2 are typically required for 300 
GeV/c 2 new quarks to form bound states, which remains 
perfectly consistent with the current bounds from the 
electroweak precision observables. As shown on Figure 
| for a new up- like fourth generation quark t' , the elec- 
troweak precision fits typically restrict the mixing ma- 
trix element \Vt>b\ — \Vtb'\ to be smaller than 10" 1 for 
m t , > 700 GeV/c 2 . If lighter, the t' total width can 
be smaller than 10 MeV [IT], thus allowing new chiral 
quarks to form bound states with nanosecond lifetimes. 




0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 



FIG. 2. Reproduction of the \ 2 fit of the oblique parameters S 
and T as given from Fig. 2 in [11] , considering a chiral fourth 
generation quark t' . The upper limits at 68% (darker blue) 
and 95% (lighter blue) CL on the t' mass are depicted as a 
function of the CKM mixing |Vt't|, freely varied over [300,800] 
GeV/c 2 and [0,0.3], respectively. The solid (dashed) curve 
denotes the upper bound m t ' (| V t >b\) under which t'i' and t'q 
(t'q) bound states can form, as defined by relation Q . 

As far as non-chiral fermions are concerned, the elec- 
troweak precision constraints provide upper bounds on 
the mixing parameters, but are not as restrictive given 
that their effects decouple in the limit of large masses 
[2"2"] . Although vector-like quarks are usually considered 
to couple dominantly to the heaviest third quark family 
(see O [23] for some exceptions), they are allowed, in 
principle, to mix with all up- (down-) type quarks |14j . 

In any case, \kqq\ < 10 -2 is used as a rule of thumb in 
most of the allowed vector-like representations 14] . Fur- 
thermore, large corrections to the CKM matrix elements 
are required to be small for SU(2)l singlets [24 , dou- 
blets and triplets [25]. In this framework, the hadronic 
production of new bound states involving sequential or 
non-chiral quarks is a conceivable possibility at the LHC. 
Assuming that fff\ is satisfied, their lifetime can be longer 



than the hadronisation timescale Aq C£) ~ 10 -23 s needed 
to build up new quarkonium resonances, open-flavour 
mesons and baryons |20) . 

If their formation is governed by the same strong in- 
teractions that are responsible for the existence of the or- 
dinary hadronic spectroscopy, the associated mass spec- 
trum can be determined from the known properties of 
low-energy hadron physics [26) . Considering a simplified 
model of QCD-like hadrons, new heavy quarks can form 
QQ bound states (hereafter defined as t/q) from induced, 
strongly attractive potentials. For instance, quarkonium 
formation can arise from Higgs boson exchange, as the 
corrections from Yukawa-type forces cannot be neglected 
for large quark masses [27[[28]. In [29], Hung detailed the 
numerical resolution of a general non-relativistic Higgs 
exchange potential, with interesting consequences on the 
associated scalar spectrum. More recently, Enkhbat et 
al. gave a preliminary discussion of the fourth genera- 
tion Yukawa bound states phenomenology at the LHC, 
considering a ultra-heavy degenerate new quark family 
[30) . Although binding energies of Yukawa origin cer- 
tainly provide a more consistent framework for a dedi- 
cated analysis, we restrict the present discussion to the 
context of Coulomb-like potential models, and assume 
that the short-distance behaviour of the quark-antiquark 
potential dominates (we refer the reader to [32] and [31] 
for previous and more exhaustive phenomenological re- 
views). In this context, the effects of the quarkonium 
wave function must be taken into account as we con- 
sider the corresponding production and decay modes. In 
perturbative QCD, the dominant contribution to i]q pro- 
duction cross-section in pp collisions reads 32 

7T 2 T dL 

o(w -> 99 -> no) = (r) Q -> gg) b~^]gg, (8) 

8a 2 (m 2 ) 

r(v Q ^99) = ^^ li \Rs(0)\ 2 , (9) 

where m v = 2niQ 7 and [T^] gg denotes the gluon lumi- 
nosity with t = ml/s. The full NLO partonic cross- 
sections for the gg, gg and qq initiated reactions are 
provided in [33], and have been updated in |34j . In- 
terestingly, they are all proportional to the square of 
the QQ radial wave function R$ at the origin, which 
might drive the heavy quarkonium production cross- 
section u(gg — > rjo) to be substantially smaller than 
a (99 ~~ ^ QQ)- Indeed, the quarkonium production rate 
typically amounts to a few percent of the pair production 
cross-section [331136) . 

Should it preferably hadronise into a quarkonium state, 
such a suppression should be taken into account when 
setting limits on a new long-lived heavy quark Q. 

Depending on the ratio of its intrinsic width Fqq to 
the associated binding energy A VQ , three distinct cases 
must be considered for a given bound state. If it is 
strongly bound such that ^qq < A w , for example 
through new confining interactions going beyond the SM, 
the corresponding bound states will appear as quarko- 



nium resonances while the constituent particles would 
not be directly observed. On the other extreme, if the 
width is much larger than the binding energy, namely, if 
Tqq > A VQ , QQ pair production could overtake t/q pro- 
duction with little evidence of any resonance behaviour 
over the continuum. The interesting and most challeng- 
ing case is certainly the intermediate regime ^qq — A,^ . 
This can occur if the new fermions hadronise through 
the strong force, while the single particle decays are sup- 
pressed. The bound states then either undergo single 
quark decays or annihilate hadronically, with specific fi- 
nal state signatures. 

We now give a brief and qualitative description of some 
of the expected r]Q final state signatures at the LHC, 
considering either a chiral or vector-like heavy quark 
Q. While similar results hold for all S- and P-wave 
quarkonium states, we limit our analysis to the case of 
a stable, neutral J PC = h pseudoscalar state r\Q. In 
particular, a ! So quarkonium state formed from down- 
type fourth generation b' quarks is known to provide a 
possible candidate if the condition ^ is satisfied [34] . 
The latter is mainly produced through the gluon fusion 
process ([£]) with a production cross-section two orders 
of magnitude larger than the J PC = 1 vector state 
[34] , If b' is the lightest fourth generation quark, a W 
bound state can decay either through q — b' family mix- 
ing with q = u, c, t, or to boson pairs. Whether the sin- 
gle quark decays would compete with the r\v hadronic 
modes depends on the heavy quark masses and mix- 
ings. If mi,' > 1 TeV/c 2 , a fourth-generation down- 
type quark dominantly decays into fermion-antifermion 
pairs if \V qb ,\ > 1(T 2 . For m b , < 1 TeV/c 2 , the rjy 
total width lies below 1 GeV/c 2 , and the quarkonium 
state decays proceed dominantly via the annihilation di- 
agrams depicted in Figure pj] [34] . Decays to gauge boson 
pairs, fermion-antifermion pairs, gauge-Higgs and Higgs 
boson pairs can occur through 7, Z and H exchange in 
the s— channel, or proceed via quark mixing in the t— 
and u— channels, if allowed. While the processes involv- 
ing charged currents are suppressed in the case of small 
quark couplings, the r\y WW mode is allowed at loop- 
level via the exchange of the SU{2)l partner t 1 , even if 
\Vtb> I — 0. Decays to WW pairs are not allowed for t\t.b 
bound states involving singlet vector-like quarks. Decays 
to Higgs boson pairs (j\q — > HH) are forbidden by CP 
conservation in the case of J PC = h pseudoscalars. 

An exhaustive study of the production and subsequent 
decays of the S— and P— wave quarkonium states of a 
heavy chiral quark is available in [52] and [37], where all 
the allowed production mechanisms and decay patterns 
have been thoroughly investigated. As an illustration, 
we compare in Figure [4] the branching ratios of the al- 
lowed decay modes for quarkonium states formed by a 
chiral fourth generation quark b' and a vector-like B_i/ S 
isosinglet [37] . The quarkonium annihilation decay rates 
depend markedly on the mass scale. 

In the chiral case, the main decay modes are ZH, gg, 
WW, ZZ , Z7 and 77 respectively. The fermionic decays 
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(d) (e) 



FIG. 3. First order diagrams for the r\Q — > p\P2 decays, where 
pi and P2 can be a gauge boson pair, a fermion-antifermion 
pair, a gauge-Higgs boson pair, or a pair of Higgs bosons. The 
strong decay mode gg — > gg is also allowed but not shown 
above. The (a), (b), and (c) diagrams denote decays via 7, Z 
or H exchange in the s— channel, while (d) and (e) correspond 
to the P1P2 decays involving quark exchange in the t— and it- 
channels. Figure from |32j . 



T]b' —> qq are also allowed if \V q b>\ 7^ 0. For pseudoscalar 
rib' masses lighter than 450 GeV/c 2 , the dominant anni- 
hilation mode proceeds via the strong interaction, i.e., 
Vb' 99- Substantial decay rates to qq fermion pairs are 
allowed in case of large Q— q couplings. For larger masses, 
the width into two gluons is decreasing, whereas the an- 
nihilation mode into ZH pairs takes over. This follows 
from the fact that the decay rates into Higgs and longi- 
tudinal gauge bosons are enhanced by the large Yukawa 
coupling of the heavy b' quark, so that the ZH branching 
ratio becomes sizeable for increasing masses. The latter 
channel then leads to a salient signature for heavy fourth 
generation bound states. Although it is suppressed below 
the percent level, the distinctive diphoton decay mode 
also provides a possible golden channel for discovery. 

While fourth generation pseudoscalar quarkonia de- 
cays allow for striking signals, TT and BB bound states 
of vector-like quark singlets can only decay to gg, 77, 7Z 
and ZZ pairs, as the t]q — > ZH mode is absent [37 . This 
follows from the fact that the ZH mode proceeds through 
the axial part of the neutral current coupling, which in 
turn is proportional to the third component of the weak 
isospin. Only 77, Z7 and ZZ would then signal the 
visible hadronic modes with branching ratios of about 
1(T 4 . TT and BB bound states of SU{2) L vector-like 
singlet quarks are thus hardly observable at the LHC as 
the gluon mode dominates. On the other hand, T and B 
vector-like doublet (triplet) partners with non-vanishing 
weak isospin can have large branching ratios to ZH via 
7, Z or H exchange in the s— channel. The rjx and rjy 




FIG. 4. Branching ratios as given in |37| for t\q quarko- 
nium, considering a chiral b' quark (top) and a down-type 
singlet vector-like B quark (bottom), with mn = 125 GeV/c 2 . 
\V t t'\ — 1 is assumed in the fourth generation case. 



quarkonia provide a noticeable exception, given that the 
exotic quarks X 5 / 3 and ^4/3 do not couple to neutral 
currents at tree-level. Since they only mix with the other 
states via charged currents, the annihilation modes in- 
cluding 7, Z and H exchange are thus forbidden for rjx 
and r/Y, which then mainly decay to gg, WW and to the 
suppressed 77 modes. 



E. Signatures for heavy stable particles : 
open-flavour mesons 

While QQ quarkonium resonances might be challeng- 
ing to observe, the possibility for "open- flavour" hadrons 
Qq (Qq) and Qqq (Qqq), with q being a light SM quark, 
also cover a broad spectrum of new heavy mesons and 
baryons to search for. Should they be produced at the 
LHC, such states can form by capturing a light quark 
partner from the vacuum. As they pass through the 
detector, they can transform into various slowly-moving 
heavy states. Assuming that they hadronise with u, d 
and s quarks after being produced, Table |V| lists the cor- 
responding stable mesons and baryons. 



The properties of new hypothetical fourth generation 
quarkonium and open-flavour mesons have been exam- 
ined previously in |38j considering large-iV expansion 
techniques. In [39], their masses and decay constants 
have been estimated from the experimental measure- 
ments of the ordinary bb (cc) hadron masses, using the 
QCD sum rules. Interestingly, both approaches indicate 
that the corresponding spectra share numerous features 
with most of the scenarios predicting stable bound states 
of heavy quarks beyond the Standard Model, depending 
on their masses, charges and angular momenta. 

If they decay within the detector, the interactions of 
such open-flavour mesons with the material occur to be 
similar to that of R- hadrons, i.e., stable hadrons com- 
posed of a supersymmetric particle and at least one SM 
quark (see [40l |4T] and references therein for related 
works). The 1/2 difference in spin with respect to stop 
hadrons has little impact on the search strategies, which 
depend mostly on the interactions between the lighter 
quarks and the detector. This is in substance the de- 
scription given by the spectator model [HI 03] , for which 
the allowed transitions of heavy and light quarks within 
a given bound state are known to be flavour and spin 
independent [44, 45j. In this context, the heavy quark Q 
can be considered as a source of kinetic energy which sole 
function is to give mass and momentum to the underly- 
ing hadron. Such new heavy bound states then behave 
as rather passive objects, consisting of a non-interacting 
heavy component, accompanied by lighter constituents 
(u and d quarks, mostly). While the partons are being 
scattered within the detector, their cross-sections vary 
with their inverse mass in perturbative QCD [6j[26]. A 
heavy Qq (Qq) bound state is then expected to suffer 
small energy losses when interacting with the detector, 
given that only its light quark half is responsible for the 
hadronic interactions with the detector. 

We list in Table IVII the various allowed 2 — !• 2 and 
2 — > 3 transitions involving meson-to-meson and meson- 
to-baryon conversions. In general, these interactions lead 
to simultaneous pion emission as the initial mesons con- 
vert into slowly-moving baryons. 

While they interact with the detector, most of the 
new states suffer multiple scattering and convert into 
baryons, allowing for Quu, Qud and Qdd states, as new 
heavy mesons are kinematically favoured to increase their 
baryon number by emitting one or more pions |26j . Due 
to the lack of these in the material, the reverse reaction 
is known to be less favoured. 

Interestingly, new mesons which convert at the be- 
ginning of the scattering chain could generate tracks 
with "dashed-lines" , signalling possible baryonic or elec- 
tric charge exchange. Indeed, Qq and Qq bound states 
traversing a medium composed of light quarks likely flip 
their electric charge, frequently interchanging their par- 
ton constituents with those of the material nuclei. Even- 
tually, they could leave the detector as heavy stable 
neutral particles and hence not be observable in muon 
detectors. The modeling of the nuclear interactions of 
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Charge 


Mesons Baryons 


= 
= 1 
= 2 
= 3 


TO, Tu, Bd, Bd, Bs, Bs Tdd, Tds, Tss, Bud, Bud, Bus, Bus, Yuu, Yuu 
Xu, Td, Ts, Bu, Yd, Ys Xdd, Xds, Xss, Tud, Tus, Buu, Bdd, Bds, Bss, Yud, Yus 
Xd, Xs, Yu Xud, Xus, Tuu, Ydd, Yds, Yss 

Xuu 



TABLE V. Qq (Qq) mesons and Qqq (Qqq) (anti-) baryons involving Q = X5/3, T2/3, -B-1/3 and YI4/3 vector-like quarks. The 
states in bold font are hadrons which yields are expected to be substantial at the LHC. For simplicity, only the neutral and 
positively charged hadrons are displayed. 



Meson-to-Meson 2 — > 2 processes: 
Charge exchange: 
Qu + p < — > Qd + n 



Elastic scattering: 
Qq + p/n^Qq + p/n 



Meson-to-Baryon 2 — > 2 processes: 
Baryon exchange: 
Qu + n — > Qud + 7r~ 
Qu + n-> Qdd + n° 
Qu + p — > Quu + n~ 
Qu + p — > Qud + 7t° 



Qd + n 
Qd + n - 
Qd + p- 
Qd + p- 



■ Qud + 7T° 

Qdd + 7T+ 
Quu + 7T° 

Qud + 7T + 



Meson-to-Meson 2 — > 3 processes: 
Inelastic scattering: 
Qu + n — > Qd + n + tt~ 
Qu + n — ¥ Qu + n + ir a 
Qu + n — > Qu + p + 7r~ 

Qli + P — > Qd + P + 7T~ 

Qu+p ^ Qu + p + n° 
Qu + p — > Qu + n + 7r + 
Qu + p — > Qd + n + 7T° 

Qd + n — > Qu + n + 7r + 
Qd + n Qd + " + t° 
Qd + n Qu + p + n° 
Qd + n — > Qd + p + 7r~ 

Qd + P -> 0" + p + vr+ 
Qd + p Qd + P + 7T 
Qd + P — > Qd + n + n + 



Meson-to-Baryon 2 — > 3 processes: 
Baryon exchange: 
Qu + n — > Qud + tt~ + it 

Qu + n — > Qdd + 7T~ + 7T + 

Qu + n -+ Qdd + 7T° + 7T° 

Qu + p — > QUU + 7T~ + 7T° 
Qu + p — > Qud + 7T~ + 7T + 
Qu + p -> Qud + 7T° + 7T° 
Qu+p^r Qdd + 7T+ + 7T° 

Qd + n — > Qud + ir~ + tv + 
Qd + n -+ Qud + 7r° + n° 
Qd + n -> Qdd + 7T° + 7r+ 

Qd + p — > QUU + 7T~ + 7T + 
Qd + p -> QUU + 7T° + 7T° 

Qd + p^ Qud + 7r° +7T+ 



TABLE VI. Meson-to-Meson and Meson-to-Baryon processes for Qu and Qd mesonic states, where Q = -X5/3, T2/3, -B-1/3, Y-4/3 
and q = u,d. Similar transitions can be obtained for the charge conjugated processes involving Qq. Assuming that the cross- 
sections for all the transitions are of the same order of magnitude, Meson-to-Baryon processes to Qud states are more likely 
than those involving Quu and Qdd in the final state. 



new heavy hadrons traveling through matter, a survey 
of which is given in 05] , actually favours scenarios 
with significant charge suppression if heavy open-flavour 
mesons have a sizeable probability to transform into neu- 
tral particles while traversing the detector. This poses a 
serious challenge for the experimental searches, since the 
most conservative scenarios assume a complete charge 
suppression where 100% of the produced hadrons become 



neutral before reaching the muon detectors. Taking a 
specific example, we consider pair-produced T 2 / 3 quarks 
hadronising into the heavy states Tq and Tq immedi- 
ately after production. For convenience, we assume that 
a majority (> 90%) of the new heavy quarks initially 
hadronise into mesons, while a smaller amount (< 10%) 
form baryons [B]. As they interact farther in the detector, 
most of the Tq mesons eventually exit the detector as Tu 
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and Td mesonic states, or as Tud baryons. Baryon-to- 
meson conversions are allowed as well, as these states can 
also annihilate back into Tu and Td mesons, yet with a 
smaller rate [17] ■ Regarding bound states of vector-like 
antiquarks, Tq mesons unlikely give rise to antibaryons, 
but can still flip their charge through the exchange of u 
and d quarks with the material. Possibly large flavour 
fractions for Tu and Td can thus be expected in the de- 
tector, with a negligible amount of Tqq antibaryons. 

After travelling through the detectors a few nuclear 
lengths away from the production vertex, the novel 
hadrons transform in roughly 100% to the positively 
charged baryon Tud, whereas one half of the T— states 
thus transform into Tu mesons, and the other half in 
Td. The corresponding hadronic flavour decompositions, 
evaluated as a function of the penetration depth in the 
detector, can be read from fig. 5 of [55] for R-hadrons 
involving stable scalar top and antitop quarks. We no- 
tice that similar results are obtained when replacing stop 
t squarks by stable quarks T of equal mass. The Tud, 
Tu (Tu) and Td (Td) states are then expected to retain 
the largest flavour composition fractions when consider- 
ing penetration depths larger than 3 meters, and are thus 
the most likely observable states, which we emphasised 
in bold font in Table EJ 

Except for Yuu and Yuu, it is interesting to notice that 
none of the states listed in Table [V] is neutral if involving 
vector-like quarks with exotic charges. The correspond- 
ing bound states might give rise to possibly large frac- 
tions of slowly moving charged particles Xud and Yud 
with Q = +2e and — e respectively, accompanied by a 
small amount of Xu, Xd, Yu and Yd mesons (and charge 
conjugates). These particles are all electrically charged as 
they go across the detector, emitting charged and neutral 
pions in the process together with losing small amounts 
of energy. Given the small expected interaction cross- 
section from the Meson-to-Baryon transitions, the pro- 
cesses shown in Table VI leave small energy deposits in 
the calorimeters, of the order of 0(1) GeV [BJ. On the 
other hand, such stable hadrons would lead to observ- 
able tracks due to the ionisation energy losses, allowing 
for signatures similar to slow-moving muons with high 
transverse momentum. Searches for stable charged par- 
ticles, if adapted to such a case, thus provide a promising 
strategy to rule out the possibility for novel exotic quarks 
with large lifetimes. Given these very specific signals, the 
aforementioned signatures can certainly be discriminated 
at the LHC. 



III. EXPERIMENTAL ASPECTS 

In this section, we review the past and current exper- 
imental searches for new long-lived quarks. Their possi- 
ble limitations are then discussed, as novel heavy quarks 
might lead to signatures involving displaced vertices, new 
hadronic bound states, or decays happening outside of 
the detector acceptance. 



A. Previous searches at Tevatron 

The searches undertaken at the Tevatron resulted in al- 
ready stringent mass bounds on long-lived heavy quarks, 
yet not without assumptions. Looking for long-lived par- 
ents of the Z boson in displaced vertices from pp col- 
lisions at s/s =1.8 TeV, CDF set limits on the cross- 
section of a fourth generation charge — l/3e quark as a 
function of its lifetime with an integrated luminosity of 
90 pb _1 [IB]. Isolated electron-positron pairs originat- 
ing from Z decays were searched for in the exclusive 
channel b' — > bZ with my > rnz- Finding no evidence 
for new long-lived particles, CDF excluded my < 148 
GeV/c 2 for cr = 1 cm at 95% CL. This limit drops to 
rn z + m b ~ 96 GeV/c 2 if ct > 22 cm or cr < 0.009 cm. 
Previously, the DO collaboration already ruled out the 
range < nty < rnz + from b' — > 67 searches 

in |49j for all proper lifetimes, while a b' quark with a 
mass lower than m^/2 was previously dismissed by the 
LEP direct searches [50) . Within a 90 pb _1 data sample 
of pp collisions recorded during 1994-95, CDF performed 
a search for low velocity massive charged stable particles 
leaving large amounts of energy in the calorimeters [51] . 

Assuming a muonlike penetration and searching for an 
anomalously high ionisation energy loss signature, the 
data was found to agree with background expectations, 
and upper limits of the order of 1 pb were derived on the 
production cross-section. Sensitive to long-lived fourth 
generation quarks scenarios, the lower bounds % > 190 
GeV/c 2 and mf > 220 GeV/c 2 have been obtained for 
q = — l/3e and q = 2/3e stable quarks respectively, with 
no observed excess over background. More recently, DO 
studied a 1.1 fb _1 data sample and looked for Z — > e~e + 
decays assumed to follow from a long-lived b' quark par- 
ent with BR(6' -> bZ) = 100% [5J]. Assuming that the 
electromagnetic showers were originating from the same 
vertex, the analysis found no hint away from the pp in- 
teraction point. DO excluded my < 190 GeV/c 2 at the 
95% confidence level for decay lengths between 3.2 mm 
and 7 m. With the same integrated luminosity, CDF 
excluded my < 268 GeV/c 2 at 95% CL, considering a 
long-lived b' quark decaying exclusively into a Z boson 
and a b jet [53] . However, it has been emphasised in 
[51] that the assumption BR(6' -> bZ) = 100% is inac- 
curate for my > 255 GeV/c 2 , given that the b' decay 
mode b' — > tW should take over if my > m t + my\r. Fur- 
thermore, the processes b' — ¥ (u, c)W were hinted to pro- 
ceed with non-negligible rates, leading to an even lower 
branching ratio. Indeed, even if the b' — » tW decay domi- 
nates, the aforementioned mass limits depend sensitively 
on the CKM mixing elements between the fourth and the 
first three generations. If non-zero couplings between the 
fourth and the lighter quarks are allowed, the CDF lower 
bound on my can be significantly affected. 

Additionally, it is known that the above CDF limits do 
not apply for long-lived heavy quarks decaying between 
roughly 1 cm and 3 m within the detector |54j . Should 
their couplings to SM quarks lie in the corresponding 
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range 4.5 x 10~ 9 — 7.8 x 10 -8 , there exists no bounds 
on the t' and b' masses in this uncovered region. For 
decay lengths larger than the detector dimensions, the 
lower bounds drop to m t ' > 220 GeV/c 2 and my > 190 
GeV/c 2 from the stable quark searches [5TJ[52]. Although 
most experiments require the quarks to decay within a 
few centimeters from the beam pipe, this is generally 
not the case for new heavy particles having very small 
mixings with the SM quarks. This issue should be treated 
carefully at the LHC. 



B. LHC 

1. Limitations of direct searches 

While dedicated searches for long-lived bound states 
of vector-like heavy quarks are still missing at the LHC, 
various searches for heavy quarks of zero lifetime have 
already been performed by the ATLAS and CMS exper- 
iments. In this section, we discuss how the results of 
searches for prompt production could be reinterpreted 
for lifetimes larger than 10~ 10 s. 

The main aspects of reinterpreting these results are the 
branching ratios to the investigated final states, which 
may be different in alternative models, and of course the 
lifetime of the heavy particle. 

Considering here a specific example, the current best 
limit on b' production, published by the CMS collabora- 
tion, excludes production cross-sections of a > 0.1 pb at 
the 95% confidence level with an integrated luminosity 
of 4.9 fb" 1 at </s = 7 TeV [55] ■ assuming a branching 
ratio BR(6' -> tW) = 100%. As we have seen in Sec- 
tion |II B| the assumption of exclusive branching ratios is 
not generally valid for new heavy quarks. Nevertheless, 
the fraction to which a hypothetical signature occurs can 
be used to recalculate the cross-section limit in the most 
naive approximation 



^truc — ^excluded ^ 



BR, 



assumed 



BR t 



(10) 



where <7 xciudcd is the excluded cross-section under the 
assumption of BR assumoc i, and BR truc is the branching 
ratio in the given model. 

While this statement is certainly trivial, the reinter- 
pretation for longer lifetimes requires more thought. In 
particular, the implicit assumption of heavy quarks de- 
caying promptly at the production vertex is valid only for 
very short lifetimes. If these particles form bound states 
and propagate certain distances before their decay, they 
can potentially escape the direct searches. In such a sce- 
nario, the published limits could be considerably weaker, 
depending on the particle's lifetimes, simply because they 
would fly too far to be detected at the primary vertex. 

For intermediate lifetimes with displaced decay vertices 
within the detector volume (cfr. the region («) discussed 
in Section II C), a recalculation of the published limits 



can be attempted. Assuming an exponential decay func- 
tion, a fraction of the decays will always happen in the 
vicinity of the beam line, so that the prompt searches will 
pick them up. Based on this fraction, one can recalcu- 
late the limits as a function of the heavy quark lifetime. 
However, in order to do this correctly, one requires the 
exact selection efficiency as a function of the displace- 
ment from the beam line. Unfortunately, such informa- 
tion has not been made available by the experiments. 
Still, it can be estimated from an educated guess under 
specific assumptions. For instance, the analysis in [55] 
applies b— jet tagging, which requires high-quality tracks 
to originate within the inner detector volume. The CMS 
innermost part, the pixel detector consists of three barrel 
layers with the innermost layer at a distance of 4.4 cm 
from the beam line. This represents the technical limi- 
tation to this analysis, so that we can assume that the 
selection efficiency vanishes at a transverse displacement 
of around 4 cm. 

Even in analyses without application of b— jet tagging, 
stringent quality cuts are usually required to select recon- 
structed detector objects originating from the primary 
vertex. One of the motivations for these requirements 
is the rejection of pileup. Jets from displaced decays 
which do not point to the production vertex are mostly 
removed by the pileup cleaning procedures. We therefore 
make the assumption that events are kept in the direct 
searches if the decay happens at less than 4 cm from 
the primary interaction vertex, and lost otherwise. The 
fraction of the lost events can be determined from the 
distance these particles travel before their decay. This 
distance d depends on the mass, the momentum and the 
proper lifetime r. It is defined in the laboratory frame as 
d = 7 • • c- t. The fraction is obtained by convoluting d 
with an exponential decay function and the momentum 
spectrum which determines 7 and /3. 

We take the momentum spectrum of pair-produced 
heavy quarks as predicted by the MadGraph [56] Monte 
Carlo generator. The assumed center of mass energy is 
yfs = 8 TeV. The momentum is needed for a precise cal- 
culation of the decay distance on an event by event basis 
because the velocity may be significantly smaller than c. 
The pt distribution for various quark masses is shown in 
Figure [5] 

Due to the very high masses these particles can be 
slowly moving. The distribution of the velocity in units 
of c is shown in Figure [6j It is visible that the average 
velocity, for instance for tuq = 1000 GeV/c 2 is about 
0.6c. 

This simulation assumes that the production mecha- 
nism and kinematics are independent of the lifetime of 
the particle. Possible limitations of that assumption are 
discussed in Section |ill B 4| Calculating the distribution 
of decay vertices in the detector frame, we then obtain 
the fraction of decays outside of the geometrical accep- 
tance of 4 cm as a function of the proper lifetime r. Table 
|VII| summarises our results for a few benchmark points 
in the parameter space. The results are displayed as a 
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FIG. 5. Transverse momentum distribution of pair-produced 
heavy quarks as predicted by MadGraph. 
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FIG. 6. Velocity distribution of pair-produced heavy quarks 
in units of /? = v/c as predicted by MadGraph. 



function of the proper lifetime in Figure [7J 



Lifetime 


10- iu s 


lO" 9 s 


10" s s 


m 
m 
m 
m 


= 300 GeV/c 2 
= 500 GeV/c 2 
= 700 GeV/c 2 
= 1000 GeV/c 2 


32% 
26% 
21% 
16% 


86.5% 
84.1% 
82.3% 
79.4% 


98.5% 
98.3% 
97.9% 
97.5% 



TABLE VII. Fraction of rejected decays of heavy quarks out- 
side of the geometrical acceptance of 4 cm around the primary 
interaction vertex for three different lifetimes and four differ- 
ent heavy quark masses. 

We can now immediately conclude that direct searches 
are only valid for lifetimes which are considerably shorter 
than about 10~ 10 s, otherwise the particles would propa- 
gate too far and be rejected by the selection criteria. We 
also stress that this result is relatively independent of the 
particle masses within our approximations. 

The published limits on the heavy quark production 



FIG. 7. Fraction of rejected decays of heavy quarks outside 
of the geometrical acceptance of 4 cm around the primary 
interaction vertex as a function of the proper lifetime r. 



cross-sections can be simply reinterpreted in the region 
l(T 8 s < t < 10- 10 s by 

Otruc = ^excluded X (1 — /) (H) 

where <J C xciuded and <7true are the published and recalcu- 
lated limits and / is the fraction of lost events given in 
Table |VII| and Figure [7J The recalculation of the mass 
limits can then be done by comparing the recalculated 
cross-sections with predicted cross-sections as a function 
of the mass. 



2. Displaced vertices 

The immediate question that arises from the previous 
results is whether lifetimes larger than 10 -10 s are cov- 
ered by dedicated searches for long-lived particles, either 
through displaced vertices or signatures of stable parti- 
cles propagating through the full detector. 

As already mentioned above, dedicated analyses for 
long-lived vector like quarks have not been done at the 
LHC, but reinterpretations of other searches may be pos- 
sible. 

Displaced topologies have been searched for at LHC 
and published, for instance in [57] and [55]. These 
searches are not directly applicable to heavy quarks for 
different reasons which will be discussed in this section. 
The sensitivity of these searches is impressive as the 
excluded cross-sections are comparable with predicted 
heavy quark production cross-sections, so it is worthwhile 
discussing this in more detail. 

The CMS analysis [57] looks for heavy, long-lived x 
bosons decaying to di-leptons x ~~ * l + l~ ■ This analysis 
is feasible thanks to the CMS track reconstruction algo- 
rithm which is able to identify very displaced tracks not 
originating from the primary interaction vertex. Tracks 
of two oppositely charged leptons are fit to a common 



13 



vertex, requiring the vertex to fulfill x 2 /n < 5 where n 
is the number of degrees of freedom. With this method 
transverse displacements up to 50 cm from the beam line 
can still be reconstructed. This analysis is able to put 
limits on production cross-sections of a x BR < 10~ 3 pb 
at 95% C.L. The b' -> tW and B tW decays, for in- 
stance, can lead to displaced di-lepton vertices as well, 
because both W bosons can decay leptonically. For a 
mass of my = 500 GeV/c 2 the predicted b'b' produc- 
tion cross-section is a = 0.310 pb [59], which is clearly in 
the vicinity of the analysis sensitivity. However, in this 
analysis, the momentum of the vertex is required to be 
parallel to the vector pointing from the primary vertex 
to the x boson decay vertex. 

This requirement is actually the limiting criterion 
which makes the result insensitive to heavy quarks. Due 
to the presence of the neutrinos and of the 6-quark jet in 
the decay cascade, the momentum of the di-lepton vertex 
will not be parallel to the vector pointing from the pri- 
mary to the secondary vertex. Consequently, di-leptonic 
decays of heavy quarks are mostly rejected in this anal- 
ysis. 

We suggest a straightforward extension of this analysis 
removing this collinearity cut. The side effect would be 
increased background levels because the purpose of this 
cut is its rejection. It is a priori not trivial to predict 
the amount of additional background from the available 
information, nor obvious whether the additional back- 
grounds can be discriminated with other methods. Given 
the high sensitivity of this analysis, as reported above, it 
should still be feasible to reach exclusions of heavy quark 
production cross-sections. 

Another component that would be necessary to derive 
limits on heavy quark production from displaced vertex 
searches is the Monte Carlo simulation of the signal pro- 
cess. The interaction of these heavy particles with the 
detector material and the efficiency of their reconstruc- 
tion and selection need to be estimated in order to facili- 
tate the calculation of limits. Difficulties may arise from 
missing implementations of certain exotic models in the 
simulation software. Existing simulations of heavy stable 
particles such as those of R-hadrons may provide good 
approximations though 26 . 

A displaced vertex analysis has also been performed 
by the ATLAS collaboration [55] , yet following a slightly 
different strategy. Displaced vertices are reconstructed 
in an inclusive way, using all displaced tracks in the 
event as potential seeds. In contrast to CMS, where only 
opposite-sign di-muon vertices have been used, the AT- 
LAS approach requires at least four tracks at the vertex 
and a vertex mass of at least 10 GeV/c 2 . Requirements 
on momentum collinearity have not been made, so heavy 
quark decays will in principle be selected. To ensure a 
good fit quality, only vertices within the fiducial barrel 
pixel detector volume are considered, which means that 
transverse displacements up to 18 cm are considered. In 
addition, one muon with pj> > 45 GeV/c is required to 
be present in the event. The muon is not required to 



originate from the displaced vertex in order to enhance 
sensitivity for models in which it arises from the proton- 
proton vertex. This analysis has been carried out in the 
context of long-lived supersymmetric particles, and limits 
on a ■ BR have been derived. Depending on particle type 
and lifetime, the limits are in the range of 0.2 to 10 pb 
using an integrated luminosity of only 33 pb^ 1 at s/s = 7 
TeV. With more data we can expect this analysis to be- 
come sensitive to heavy quark decays as well. However, 
no attempt has been made to put this into the context 
of heavy quark searches. So the obvious suggestion is to 
extend the theoretical models covered by this analysis to 
heavy quarks as well. Also in this case the Monte Carlo 
simulation of the signal process is the central ingredient 
for an extension of the coverage to heavy quarks. 



3. Heavy Stable Charged Particles 

For very long-lived scenarios, the new heavy states do 
not decay, but travel through the full detector. In the 
following we review recent results by CMS [60 along with 
an assessment of their relevance for HSCPs. 

The main identification criteria for HSCPs are high 
momentum, high ionization and long time of flight 
(TOF). The energy loss dE/dx along the track is cal- 
culated from the charge deposits in the silicon tracker, 
while the TOF is obtained from the arrival time in the 
muon system. These quantities are uncorrelated for SM 
particles. This noncorrelation is used to estimate the 
background yields in signal-depleted regions. 

Models of charge suppression due to interaction of 
HSCPs with the detector material are considered as well. 
If charged particles become neutral while propagating 
through the detector they do not reach the muon cham- 
bers. The TOF criterion cannot be used in this case. 
Therefore, the results have also been estimated without 
the TOF requirement at the cost of weaker exclusion lim- 
its. 

The obtained limits for a given mass can be very dif- 
ferent depending on the assumed particle type. This is 
due to the predicted kinematics of the particle and its 
expected detector signal. For instance, the limit on pair 
production of a scalar top quark with a mass between 300 
and 800 GeV/c 2 is a < 3 fb for an integrated luminosity 
of 5 fb" 1 at y/s = 7 TeV. 

Limits for long-lived heavy quarks have not been con- 
sidered explicitly, unfortunately. However, as an approxi- 
mation, we make the assumption that a long-lived bound 
state of a stop would behave similar or equal to a bound 
state of a heavy up-type quark [25]. This concerns both 
the interaction with the detector material and the pro- 
duction mechanism which is assumed to proceed via the 
strong interaction. It is then possible to use the stop 
exclusion limits directly for stable heavy quarks. Limita- 
tions in that assumption are discussed in Section |HI B 4| 

The predicted production cross-sections of heavy 



quarks are quite large for low masses (Section II A I 
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Even with a large systematic safety margin, these cross- 
sections can be considered to be excluded by the HSCP 
searches. Extremely heavy quark masses of the order of 
rriQ w 1 TeV/c 2 are not yet excluded though, neither 
in direct searches for prompt decays nor in long-lived 
searches. 

To make a quantitative statement about the excluded 
lifetimes by the HSCP searches, we repeat our simulation 
from Section III B 1 To be detected by these analyses, 



the particles have to traverse the full tracking devices, 
and the muon chambers for the combined tracker+TOF 
analysis. The CMS tracker has a length of about 6 m and 
a diameter of 2.8 m [ST]. We can calculate the fraction of 
decays inside the tracker volume and assume that these 
decays will not be reconstructed in the HSCP analyses 
because of missing hits in the outermost tracker layers. 

Tables VIII and Figure [8] summarise our estimations 
for a few benchmark points and as a function of the 
proper lifetime. From these results we see that heavy 
quark masses reaching itiq — 600 GeV/c 2 can be ex- 
cluded at 95% C.L. for lifetimes longer than 2 • 10" 9 s. 
The transition between the two extrema of being fully 
efficient and losing all decays happens within two orders 
of magnitudes of the lifetime between 10~ 9 s and 10~ 7 s. 



Lifetime 2 


10" 9 s 


10" s s 


10"' s 


m = 300 GeV/c 2 
m = 500 GeV/c 2 
m = 700 GeV/c 2 
m = 1000 GeV/c 2 


90% 
94% 
96% 
98% 


47% 
53% 
58% 
64% 


7% 
8% 
10% 
11% 



TABLE VIII. Fraction of rejected decays of heavy quarks 
within the geometrical acceptance of the CMS tracking de- 
tector. Decays within the tracker are rejected because the 
particle's trajectory does not reach the outermost tracker lay- 
ers. 
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FIG. 8. Fraction of rejected decays of heavy quarks within 
the volume of the CMS tracking detector as a function of the 
proper lifetime r. 



4- Discussion and alternative search topologies 

We saw from the previous sections that direct searches 
for prompt decays, which exclude masses up to tuq = 600 
GeV/c 2 , are valid for lifetimes shorter than 10" 10 s. 

On the other hand, searches for HSCP can access life- 
times larger than 2 • 10" 9 s for the same mass region. 
This means that the small window in between these two 
extrema, corresponding to the region [in) discussed in 
Section II C[ cannot be excluded yet. In principle, this in- 
visible region could be probed by displaced vertices anal- 
yses, but the published analyses appear to not be sensi- 
tive enough yet for the heavy quark case. As mentioned 



in Section III B 2 displacements beyond 50 cm cannot be 
detected so that these searches are not able to cover the 
full tracker volume. 

To understand the potential of combining all three 
methods, namely, prompt searches, displaced vertices 
and HSCP searches, we repeat the simulation from above 
and calculate the fraction of decays that occur in the 
invisible region. This is the region with displacements 
larger than 50 cm, but decaying within the tracker vol- 
ume. The result is shown in Figure [9| We see that in the 
worst case scenario, 50% of the decays are lost in these 
analyses for a lifetime of r = 3 • 10" 9 s. This represents 
a moderate decrease of the selection efficiency for a nar- 
row window, so that we can expect a full coverage of the 
lifetime spectrum in the future. 
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FIG. 9. Fraction of decays of heavy quarks between the ac- 
ceptance of the displaced vertex analysis and HSCP analyses. 

To conclude, we now discuss some of the assumptions 
made in this section, and their possible limitations. In 
particular, we considered a general framework in which 
the new heavy quarks were pair-produced, with kinemat- 
ics independent of the lifetime of the particle. Our mo- 
tivation was that the current searches for long-lived R- 
hadrons formed from stable stops generally assume that 
strong pair production dominates. In supersymmetric 
models, however, the situation changes drastically when 
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the squark masses become large. As pointed out in [52] . 
the gluon fusion mechanism gets suppressed compared to 
the squark-squark and squark-antisquark pair production 
processes for large masses due to the decreasing proba- 
bility for the gluons to carry large fractions of the pro- 
ton momentum. For scenarios where either the squark 
masses or the mass of the gluino propagator becomes 
close to 1 TeV, Higgs and squark weak production start 
to dominate. The production channels then become more 
and more model-dependent for increasing masses, argu- 
ing for dedicated MC simulations above the TeV scale. 
Because of the lower production cross-section, models in- 
volving such new heavy quarks could be very challenging 
to detect at the LHC, depending on how well the signal 
signatures can be discriminated. 

Another important limitation that concerns the future 
HSCP searches for even higher masses is the required 
timing to detect such slowly moving particles. While the 
LHC design bunch spacing is 25 ns, the 2011 and 2012 
LHC runs have been using a bunch spacing of 50 ns. A 
recorded event has to be contained in a single bunch, 
even tighter constraints may be imposed by the readout 
of the tracking devices. To be visible in HSCP searches, a 
particle must traverse the full tracker and muon systems 
in order to trigger the event. Therefore, this is required 
to happen within a time frame of O(10~ 8 s). 

As discussed in the previous sections and visible in Fig- 
ure [6] the average velocity reduces with the mass. Since 
particles with a mass mq = 1000 GeV/c 2 have an av- 
erage velocity of 0.6c, this means that they propagate 
a distance of 3 meters on average within 10~ 8 s, just 
enough to pass through the tracker volumes. 

This estimation shows that the HSCP searches have 
an intrinsic limitation as the excluded masses increase 
beyond the TeV scale. Stable quarks with larger masses 
might evade the current searches as such slow states can- 
not be detected. There are not many options to overcome 
this problem. One possibility for mitigation is the search 
for short tracks. An event triggered in the muon sys- 
tem by a slowly moving particle may still contain a few 
tracker hits in the outermost tracker layers. This would 
require new track reconstruction algorithms searching for 
hits of highly ionizing particles pointing towards the trig- 
gering signature. One can measure the time of flight from 
the outermost tracker layers to the triggering muon sta- 
tions which should be large and facilitates background 
suppression. 

In case the heavy quark decays within the tracker vol- 
ume, the signature to be searched for may be one or more 
displaced jets, which trigger the event, and originate from 
the heavy quark decay. In addition a short track with a 
few hits and large dE/dx pointing to the decay vertex 
could be required. This signature may have relatively 
low backgrounds as the jets should have a large trans- 
verse momentum. The tracks within these jets arc not re- 
quired to arise from the primary interaction vertex which 
makes track reconstruction difficult. To circumvent this, 
one could perform a regional search for short tracklets 



pointing towards the energy deposit in the calorimeter. 

Besides the possible differences between the kinematics 
of long-lived quark scenarios, a final remark to this dis- 
cussion concerns their expected production rates. While 
we assumed above that the heavy quark pair production 
cross-sections were of the same magnitude as the rate of 
particles hadronising into bound states, we stressed in 
Section |II D| that this is not the case in general. As a 
matter of fact, smaller cross-sections are expected due to 
the possibly large suppression arising from the quarko- 
nium wave function, which also modifies the associated 
decay phenomenology. In particular, t]q production oc- 
curs to be smaller by at least one order of magnitude with 
respect to QQ pair production when accounting for the 
attractive potential between the two heavy quarks [34j . 
Precise calculations would be needed to refine our anal- 
ysis and the above discussion might need to be adjusted 
in case of strongly bound states. 



IV. SUMMARY AND CONCLUSIONS 

The current searches for heavy quarks at the LHC 
mostly ignore the option of long lifetimes. Although 
a new chiral family of fermions is now strongly dis- 
favoured from the recent Higgs search results, such a 
scenario remains of prime importance when considering 
vector-like quarks beyond the Standard Model. As non- 
chiral fermions decouple in the limit of vanishing mixing 
with the lighter generations, they can be long-lived and 
form bound states. Dedicated searches for novel heavy 
QQ quarkonia, Qqq baryons and Qq (Qq) open-flavour 
mesons might provide promising strategies for future in- 
vestigations. 

As we have seen, searches for prompt decays are still 
valid for long lifetimes to a certain extent, but the as- 
sumptions about branching ratios and the potential loss 
of efficiency need to be considered. We showed that these 
analyses start to become insensitive for lifetimes longer 
than 10~ 10 s, corresponding to Q — q quark couplings be- 
low the 10~ 9 level. For lifetimes beyond 10~ 10 s, searches 
for displaced vertices in the tracker volume may play an 
important role. Still, the published analyses appear to 
not be sensitive to long-lived heavy quarks, either due to 
the selection criteria that reject them or due to a lack 
of statistics. We proposed some improvements to these 
analyses to extend their reach to long-lived heavy quarks 
as well. 

Moreover, our reinterpretation of the HSCP searches 
indicates that quark masses of 600 GeV/c 2 can be ex- 
cluded at 95% CL for lifetimes longer than 2 • 10~ 9 s. 
This result however assumes that the interactions of the 
heavy quarks with the detector material are similar to 
that of stop particles. As discussed in this review, there 
are potential restrictions to such an assumption, and the 
forthcoming searches should be interpreted in a model- 
specific context. Additionally, it is usually considered 
that the bound states production cross-sections are com- 
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parable to that of pair-produced heavy quarks. This is 
not necessarily the case for the reasons we discussed. 

For long-lived particles above the TeV scale, further 
limitations include the fact that slowly moving particles 
do not fit into the nanosecond time frame required to 
read out the detector within one LHC bunch crossing. 
We suggested extensions to the existing searches, such as 
to include short tracks or alternative topologies to extend 
their reach to extremely high masses. 
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